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A new dynamic line, the Frenkel line, has recently been proposed to separate the supercritical 
state into rigid-liquid and non-rigid gas-like fluid. The location of Frenkel line on the phase diagram 
is unknown for real fluids. Here, we map the Frenkel line for three important systems: CO2, H2O 
and CH4. This provides an important demarcation on the phase diagram of these systems, the 
demarcation that separates two distinct physical states with liquid-like and gas-like properties. We 
find that the Frenkel line can have similar trend as the melting line above the critical pressure. 
Moreover, we discuss the relationship between unexplained solubility maxima and Frenkel line, and 
propose that the Frenkel line corresponds to the optimal conditions for solubility. 


Recently, there has been a dramatical increase of using 
supercritical fluids in extraction and purification appli¬ 
cations, including in food, nuclear waste, petrochemical 
and pharmaceutical industries mi- Supercritical fluids 
attract significant attention due to their extremely good 
dissolving power and “tunable” properties. The solubil¬ 
ity of supercritical fluids depend on density and diffusiv- 
ity. Supercritical fluids combine the best of both worlds: 
high density of liquids and large diffusion constants of 
gases. Moreover, both of those properties can be tuned 
over a wide range pressure and temperature above the 
critical point, optimizing their dissolving ability. 

Carbon dioxide, water and methane are three most 
commonly used supercritical fluids. In particular, H 2 O 
and CO 2 , are both abundant, non flammable and non 
toxic. They are also “non-polar” and “polar” solvent, re¬ 
spectively, so they can dissolve “polar” and “non-polar” 
solutes, respectively. The critical temperature (Tc) of 
CO 2 is at 304 K, which is near the room temperature, 
and the critical pressure (Pc) is 74 bar, which is also ac¬ 
cessible. Additionally, CO 2 can be used with co-solvents 
to modify it into ’’polar” solvent. 

The solubility of variety of solutes have been measured 
in supercritical CO 2 near the Tc as a function of pres¬ 
sure [T]. Interestingly, the experiment show intriguing 
solubility maxima above critical temperature: solubility 
first substantially increase with pressure, followed by its 
decrease at higher pressure [SHig. This effect is not cur¬ 
rently understood theoretically. Understanding it would 
lead to more efficiently use of supercritical fluids. More 
generally, it is often acknowledged that wider deployment 
of supercritical fluids and optimizing their use would ben¬ 
efit from a theoretical guidance mm- 

Until recently, supercritical state was believed to be 
physically homogeneous, which means that moving along 
any path on a pressure and temperature above the crit¬ 
ical point does not involve marked or distinct changes. 
The Frenkel line has recently been proposed, which sepa¬ 
rates two dynamically distinct states: the gas-like regime 
where particle only have diffusive motion and the liquid¬ 
like regime where particle combine both solid-like quasi¬ 
harmonic vibrational motion and gas-like diffusive mo¬ 


tion [I3HIS]. This transition take place when liquid 
relaxation time t approaches Debye vibration period, 
tq. Liquid relaxation time is defined in a usual way 
as the average time between two consecutive diffusion 
events (molecular rearrangements between two quasi¬ 
equilibrium positions) in the liquid at one point in space 
[16]. When r « td, the system loses the ability to sup¬ 
port shear modes at all available frequencies, up to Debye 
frequency, and retains gas-like diffusive dynamics only. 
The Frenkel line starts from 0.7-0 .8 Tc at Pc and ex¬ 
tends to arbitrarily high pressure and temperature on 
the phase diagram m- There are many ways to locate 
the Frenkel line on the phase diagram, yet velocity au¬ 
tocorrelation function (VAF) provides a convenient and 
mathematically meaningful criterion: the disappearance 
of oscillations and minima of the VAF correspond to pres¬ 
sure and temperature of the Frenkel line m- 

In this paper, we map the Frenkel line on the phase 
diagram using Molecular Dynamic (MD) simulation by 
calculating VAF. We study the location of Frenkel line 
for CO 2 , H 2 O and CH 4 on phase diagram, particularly 
addressing the slop of the Frenkel line in relation to the 
melting line. We subsequently compare the Frenkel line 
with the solubility maximum from experiment [MZl IS] 
and discuss why the Frenkel line is related to the solubil¬ 
ity maxima. 

We use DL_POLY MD simulation package [18] . 
and have simulated 4576 CO 2 molecules, 3375 H 2 O 
molecules and 3375 CH 4 molecules using constant- 
pressure-temperature ensemble. The intermolecular po¬ 
tential for CO 2 is the rigid-body non-polarizable poten¬ 
tial based on a quantum chemistry calculation, with the 
partial charges derived using the Distributed Multipole 
Analysis method [19]. The intermolecular potential of 
H 2 O is T1P4P/2005, which can describe intermolecular 
force very well [SS] • The intermolecular potential of CH 4 
is taken from Refs m- This potential also shows good 
accuracy in the supercritical state. We used cut-off as 12 
Afor potential, and the Smoothed Particle Mesh Ewald 
for long-range forces. We first equilibrate the system dur¬ 
ing 10 ps, and ensure the equilibration at given (P, T) 
conditions during the subsequent 40 ps. We collect and 
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analyse the result during following 50 ps. In the range 
of our simulations, the difference between MD and ex¬ 
perimental density from the NIST data base [ 22 ], is less 
than 5%. Our pressure range extends to several GPa and 
includes the pressure used in industrial applications. 

It is well known that VAF is a monotonically decay¬ 
ing function in the gas state, whereas it shows damped 
oscillation in the liquid and solid state. VAF is defined 
as: 

Z{t) =< v{0) ■ v{t) > (I) 

In the previous work m, it was shown that the min¬ 
imum of VAF would disappear when the system crosses 
the Frenkel line in the supercritical state. In Fig. 1, 
we show representative VAFs for CO 2 , H 2 O and CH 4 at 
900 bar. We can clearly see that as the temperature in¬ 
creases, the minimum becomes more shallow and finally 
disappears, which corresponds to the loss of oscillatory 
component of molecular motion and gives iP,T) for the 
Frenkel line. We note that the longitudinal mode persists 
above the Frenkel line, albeit starts disappearing with 
temperature, starting with the shortest wavelengths jSS] . 

In Fig. 2, we map the Frenkel line for CO 2 , H 2 O and 
CH 4 using the VAF criterion (disappearance of the first 
minimum). For technologically important CO 2 and H 2 O, 
we show the Frenkel line in both (pressure, temperature) 
and (density, temperature) coordinates. We also show 
the melting line [23ll26] on the phase diagram. 

We observe that the Frenkel line for all three fluids 
starts from 0.7~ 0.8 at Pc, which agrees with our pre¬ 
vious result on Lennard-Jones fluids m- Notably, the 
Frenkel line does not need to start from the critical point 
because fundamentally it is related to critical phenom¬ 
ena, and exists in systems such as the soft-sphere system 
where the boiling line and the critical point are absent 
altogether Hzj. 

We now discuss a relationship between the Frenkel line 
and the melting line, the relationship that can serve as 
a useful guide to map the Frenkel line on the phase dia¬ 
gram for any system. As discussed above, the Frenkel line 
starts slightly below the critical point. At higher pres¬ 
sure, we can predict that the Frenkel line is parallel to the 
melting line in the log-log plot on the basis of the well- 
known scaling argument. Indeed, starting from high GPa 
pressures, the intermolecular interaction is reduced to its 
repulsive part only, whereas the cohesive attracting part 
no longer affects interactions (at low pressure, the paral¬ 
lelism between the two lines holds only approximately be¬ 
cause the interactions are not well approximated by sim¬ 
ple repulsive laws, see below). In a sufficiently wide pres¬ 
sure range, the repulsive part can be well approximated 
by several empirical interatomic potentials such as the 
Buckingham-type functions or Lennard-Jones potentials 
with inverse power-law leading terms at short distances 
[/ oc ^ . For the inverse power law, a well-known scaling 
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FIG. 1: (Color online)Velocity auto-correlation functions for 
CO 2 (a), H 2 O (b) and CH 4 (C) showing the crossover at the 
Frenkel line at 900 bar. 
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FIG. 2: (Color online)The Frenkel line for CO 2 , H 2 O on the pressure-temperature (left) and density-temperature (right) 
phase diagram. The Frenkel line for CH 4 is shown in the pressure-temperature diagram. The solubility maximum of different 
solutes in supercritical CO 2 are shown in graph (a). The open circle are the solubility of /3-carotene [9]; the squares are 
l,4-bis-(octadecylamino)-9,10-anthraquinone [^; triangles are l,4-bis-(n-alkylamino)-9,10-anthraquinone [B]; the diamonds are 
biphenyl [7]; the pentagon are adamantane m-, the cross are l,4-bis-(hexadecylamino)-9,10-anthraquinone [11] 
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of pressure and temperature exists: the system proper¬ 
ties depend only on the combination of TP^ , where 7 is 
uniquely related to n. Consequently, TP^ = const, on 
all (P, T) lines where the dynamics of particles changes 
qualitatively, as it does on both the melting line and the 
Frenkel line. This implies that the Frenkel and melting 
lines are parallel to each other in the double-logarithmic 
plot, the insight that we have recently used to construct 
the Frenkel line for molecular hydrogen . 

Although our simulations were in the practically use¬ 
ful range of pressure and did not extend to high enough 
pressure to meet the condition above, we observe that 
the Frenkel line has similar trend as the melting line: for 
CO 2 , the slopes of the Frenkel line and the melting line 
both starts to increase around 1,000 bar. For H 2 O, both 
lines are flat below 1,000 bar, but their slopes start In¬ 
creasing at higher pressure. We also observe similar slope 
increase for methane simultaneously around 1,200 bar. 

The speed of sound Is one of Important properties 
shows qualitatively changes In the supercritical state. 
Notably, the speed of sound decreases with temperature 
below the Frenkel line, as in liquids and solids, but in¬ 
creases with temperature sufficiently above the line, as 
in gases [HHIS]. We note that the minima of the speed 
of sound are not absolute in the sense that their posi¬ 
tions depend on the path on the diagram (the position of 
the minimum along isobars, isochors and isotherms can 
be different). The scaling argument above implies that 
if the minima of the speed of sound correspond to the 
qualitative change of particle dynamics, the line of these 
minima should be approximately parallel to the Frenkel 
line. In Fig.2, we show pressure and temperature that 
correspond to the minimum of the speed of sound as de¬ 
duced from the NIST database [25]. We observe that the 
line of speed of sound minimum is approximately parallel 
to the Frenkel line at high pressure as predicted. 

To discuss the relationship between the solubility and 
the dynamic property of supercritical fluid, we show 
isothermal solubility maxima of different solutes in CO 2 
[SHII], on the phase diagram (Fig. 2(a)). Importantly, 
we observe the points of solubility maxima are close 
to the Frenkel line. The solubility of maxima of sev¬ 
eral solutes, such as /3-carotene, l,4-bis-(n-alkylamlno)- 
9,10-anthraquinone and l,4-bis-(hexadecylamino)-9,10- 
anthraqulnone coincide with the Frenkel line. 

We now explain the proximity of the solubility max¬ 
ima and Frenkle line, as follows. Let us fix a temper¬ 
ature above the critical point to the left of the Frenkel 
line and increase the pressure (moving horizontally to the 
right in Fig. 2). Pressure has two competing effects on 
diffusion. On one hand, it increases density and hence 
the contact area and cleaning (dissolving) efficiency. On 
the other hand, the density increase results In decreas¬ 
ing the diffusion constant. Indeed at the Frenkel line, 
where the molecular dynamics acquires the oscillatory 
component, molecular rearrangements become markedly 


less frequent. In contrast to the gas-like dynamics above 
the line where the oscillatory component of motion is 
absent. Therefore, at the Frenkel line, the supercritical 
fluid has maximal density possible at which the diffusion 
Is still in the fast gas-like regime and not in the slow 
liquld-llke regime. The optimum combination of these 
two properties gives solubility maxima. 

The data for solubility maximum for H 2 O and CH 4 
are not available. From the Fig. 2 (b) and (c), we see 
the reason why it is difficult to perform these In experi¬ 
ments. In the case of water, the pressure of Frenkel line 
Is ^9,000 bar at Tc, which is ^40 times of Pc of H 2 O 
(220.64 bar). Although, the pressure of CH 4 is not too 
high (~600 bar) at Tc, the pressure increase to ^4,000 
bar at room temperature. In both cases, their pressure 
are much higher than Pc. Compared with them, CO 2 is 
located ~ 1,000 bar at Pc, which is relatively more afford¬ 
able in experiment. We propose the Frenkel line serve as 
a predictive tool to locate the solubility maxima on the 
phase diagram. This provides a useful guide for future 
experiments. 

We note that the increase of pressure along the Frenkel 
line results in several other Interesting and potentially 
Important effects such as the increase of fluid density 
and diffusion constant as well as the appearance of the 
viscosity minimum [iMniiiii- In addition, surface ten¬ 
sion tends to zero around and above the critical point so 
that that the problem of wetting is avoided. Accordingly, 
these conditions may favour other important properties 
of supercritical fluids: for example, the speed of chemi¬ 
cal reactions may have a maximum close to the Frenkel 
line. In this and other cases, supercritical technology will 
further benefit from theoretical guidance and receive an 
impetus for using the supercritical fluids in the hitherto 
unknown range of 1-10 kbars. 

In summary, we mapped the Frenkel line for three im¬ 
portant system: CO 2 , II 2 O and CH 4 . This provides an 
important demarcation on the phase diagram of these 
systems, the demarcation that separates two distinct 
physical states with llquid-like and gas-llke properties. 
We proposed that the Frenkel line can serve as a impor¬ 
tant guide to estimate the location of solubility maxima, 
so that the cleaning and dissolving abilities of the super¬ 
critical fluids are optimized. 
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Brazhkin Is grateful to RSF (14-22-00093) for the finan¬ 
cial support. 
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